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Abstract 

The model-independent method based on the triangle ideology is implemented to 
extract the CKM-matrix angle 7 in the decays of doubly heavy baryons containing the 
charmed and beauty quarks. We analyze a color structure of diagrams and conditions 
for reconstructing two reference-triangles by tagging the flavor and CP eigenstates of 
mesons in the fixed exclnsive channels. The characteristic branching ratios 
are evaluated in the framework of a potential model setting a parametric dependence 
on the hadronic matrix elements for the decay rates. 

1 Introduction 

The current success in the experimental study of decays with the CP-violation in the gold- 
plated mode of neutral 5-meson by the BaBar and Belle collaborations P] allows one to 
extract the CKM-matrix angle (3 in the unitarity triangle by the model-independent method. 
The intensive efforts are intended in the physical programs on the B and Bg mesons at the 
hadron colliders both the active Tevatron [2j and prospective LHC. Due to the relatively high 
cross-sections the doubly heavy hadrons such as the B^ meson and baryons Pfec and Sec, 
Dftc would be copiously produced at such the machines mil In addition to the indirect or 
model-dependent measurements of unitarity triangle [S], there is an intriguing opportunity 
to extract the angle 7 in the model-independent way using the strategy of reference triangles 
m in the decays of doubly heavy hadrons. This ideology for the study of CP-violation in B^ 
decays was originally offered by M.Masetti [7j and investigated by R.Fleischer and D.Wyler [H]. 
In this letter we extend the method to study the decays of doubly heavy baryons containing 
the charmed and beauty quarks. 

To begin, we mention the necessary conditions for extracting the CP-violation effects in 
the model-independent way. 

1. Interference. The measured quantities have to involve the amplitudes including both 
the CP-odd and CP-even phases. 
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2. Exclusive channels. The hadronic hnal state has to be hxed in order to isolate the 
dehnite matrix elements of CKM matrix, which can exclude the interference of two CP- 
odd phases with indehnite CP-even phases due to strong interactions at both levels of 
the quark structure and the interactions in the hnal state. 

3. Oscillations. The dehnite involvement of the CP-even phase is ensured by the oscillations 
taking place in the systems of neutral B or D mesons, wherein the CP-breaking ehects 
can be systematically implemented. 

4. Tagging. Once the oscillations are involved, the tagging of both the havor and CP 
eigenstates is necessary for the complete procedure. 

The gold-plated modes in the decays of neutral B mesons involve the oscillations of mesons 
themselves and, hence, they require the time-dependent measurements. In contrast, the 
decays of doubly heavy hadrons such as the Be meson and baryons with the neutral 
or meson in the hnal state do not require the time-dependent measurements. The triangle 
ideology is based on the direct determination of absolute values for the set of six decays: the 
decays of baryon in the tagged meson, the tagged meson, the tagged CP-even state, 
and those of the anti-baryon. To illustrate, let us consider the decays of 

and El ^ 


The corresponding diagrams with the decay of 6-quark are shown in Fig. ^ We stress that 
two diagrams of the baryon decay to has the additional negative sign caused by the Pauli 
interference of two charmed quarks, while the color factors is analyzed in the next section. 



Figure 1: The diagrams of 6-quark decay contributing to the weak transitions El D^E^ 
and El ^ D^El 


The exclusive modes make the penguin terms to be excluded, since the penguins add an 
even number of charmed quarks, i.e. two or zero, while the hnal state contains two charmed 
quarks including one from the 6 decay and one from the initial state. By the same reason the 
diagrams with the weak scattering of two constituents, i.e. the charmed and beauty quarks 
in the El baryon, are also excluded for the given hnal state (see Fig. |2I). 

The weak scattering of 6 quark oh the charmed quark in the initial state can contribute 
in the next order in as shown in Fig. |21 Nevertheless, we see that such the diagrams have 
the same weak-interaction structure as at the tree level. Therefore, they do not break the 
symmetries under consideration. The magnitude of tts-correction to the absolute values of 
corresponding decay widths is discussed in Section |21 
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Figure 2: The penguins and weak scattering diagrams. 



Figure 3: The diagrams for the weak scattering of b and c quarks contributing to the transition 

-0 , 7 =) 0'::0 
^bc ^ ^c- 


Thus, the CP-odd phases of decays under consideration are determined by the tree-level 
diagrams shown in Fig. ^ Therefore, we can write down the amplitudes in the following 
form: 




:0 

■'6c 




:0 

■'6c 


^^^D^)=AD = Vsy:,-M 


^D, 


( 1 ) 


where Aijj D denote the CP-even factors depending on the dynamics of strong interactions. 
Using the dehnition of angle 7 


def 

7 = —arg 




for the CP-conjugated channels^ we hnd 




:0 

“'6c 


; 0 ^ 0 ) _ A{1 


:0 

'be 


= A 


D- 


( 2 ) 


We see that the corresponding widths for the decays to the flavor tagged modes coincide 
with the CP-conjugated ones. However, the story can be continued by using the dehnition of 
CP-eigenstates for the oscillating <-»• Z)° system^, 

Di,2 = T(D0±50), 

so that we straightforwardly get 

= V2Adi = Af) + Ad, (3) 

_ 72.4(5°, ^ E«D ,) g V2A'=‘f = e^‘-'At, + Ao. (4) 

^For the sake of simplicity we put the overall phase of arg VebV*^ = 0, which corresponds to fixing the 
representation of the CKM matrix, e.g. by the Wolfenstein form [H|. 

^The suppressed effects of CP-violation in the oscillations of neutral D mesons are irrelevant here, and we 
can safely neglect them. 


3 























The complex numbers entering (H and dH) establish two triangles with the definite angle 27 
between the vertex positions as shown in Fig. Thus, due to the unitarity, the measurement 



Figure 4: The reference-triangles. 


of four absolute values 


\An\ = \AiEl^ElD% 


( 5 ) 


can constructively reproduce the angle 7 in the model-independent way. 

The above triangle-ideology can be implemented for the analogous decays to the excited 
states of charmed hyperons in the final state. 

The residual theoretical challenge is to evaluate the characteristic widths or branching 
fractions. We address this problem and analyze the color structure of amplitudes. So, we find 
that the matrix elements under consideration have the same magnitude of color suppression 
A ~ 0(1/\/Nc), while the ratio of relevant CKM-matrix elements. 



~ 0 ( 1 ) 


with respect to the small parameter of Cabibbo angle, A = sin^c;, which one can easily find 
in the Wolfenstein parametrization 



( Vud 

Vus 

Vub \ 


/ 1-^A2 A AX%p-irj)\ 

VcKM = 

Vcd 

Vcs 

Vcb 

= 

-A 1 - 1A2 AA^ 


\ Vtd 

Vts 

Vtb / 


V AA^(1 — p — ir]) —AX^ 1 / 


Thus, we expect that the sides of the reference-triangles are of the same order of magnitude, 
which makes the method to be a realistic way for extracting the angle 7 . 

In Section 121 we classify the diagrams for the decays of doubly heavy baryons and 12°^ 
by the color and weak-interaction structures. Section Elis devoted to the numerical estimates 
in the framework of a potential model. The results are summarized in Conclusion. 
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2 Color structures 


Let us remind a general framework of the 1/iVc-expansion. Its physical meaning at iVc —>■ cxd 
quite reasonably implies that the quarks are bound by the gluon string, which is not broken 
by the quark-antiquark pair creation. So, the excitations of the ground states are quasi-stable 
under the strong interactions producing the decays described by the 1/iVc-suppressed terms. 
In this method we have got the following scaling rules of color structures: 


1. The meson wavefunction 


2. The baryon wavefunction 


3. The coupling constant 


4. The Casimir operators 



4/ ^ ~ 


1 


a. ~ 


Nr_ 


N‘^ -1 

Ca = N,, Cf = ^^^0{N,). 

5. The Fierz relation for the generators of SU(iVc) group in the fundamental representation 

j.Ak _ ^ I A i j.Ak 

r m — 0 rnO j I ml j- 


6 . The baryon structure constant 


Cb = 


iVc + 1 

2iV, 


~ 0 ( 1 )- 


It gives the color factor emerging in the connection of two quark lines entering the 
baryon by the gluon line (see Fig. 0. 



Figure 5: The connection of baryon structure constant. 


5 











Next, the non-leptonic weak Lagrangian has a form typically given by the following term 

mo: 

w., = ^ v^{i‘r,c,) v;,(u'‘n,) c± {s‘, s\ ± 6\ s^) + ... (6) 

where = 7^(1 — 75 ), and the Wilson coefficients 

Ci ~ 0(1) 


in the 1 /W-expansion. 

Then, we can proceed with the analysis of decays under consideration. 



All three diagrams shown in Fig. ^have the same order in 1/Nc, i.e. 


A.I ~ A.2 ~ ^3 ~ 


1 

7k' 


More definitely we get the color factors 

= = ^ = (ai - as) 

where 

[0+(W + l) + 0_(W-l)], 
[0+(W + l)-0_(iV,-l)]. 


CLi 

C 12 


2 W 

1 


(7) 

( 8 ) 
(9) 


Thus, we have to calculate the three diagrams given above in the leading order in l/W- 


2.2 




be 


The decay modes 


oo 

“fee 


and 


are described by the diagrams shown in Fig. IHl similar to those of Fig. ^ The only difference 
is the replacement of d quark by the strange one, that should be taken into account by the 
anti-symmetrization of wavefunction in the final state, i.e. the baryon structure of which 
results in the vector-spin state of the doubly strange diquark. Then, the appropriate color 
factors are given in 0 . 

The diagrams for the decay modes 


and ^ D 


O'-O 
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Figure 6: The diagrams of 6-quark decay contributing to the weak transitions 
and S», ^ D^nl 


can be obtained from those of Fig. IHlby the replacement of weak currents u —> s by u —>■ d 
and c —>■ s by c —>■ d. The triangle ideology is effective in the case under consideration too. 
The color factors are given by m again. However, the amplitude of is suppressed 

by the CKM-matrix factor of 

vcv:, 

which implies that the corresponding side of reference-triangle will be much less than another. 
In practice, the relatively large branching of decay to the meson should be measured with 
extremely high accuracy in order to make a sense in the reconstruction of the triangle with 
the relatively small side determined by the branching of decay to meson. 

However, with an expected event rate discussed below, in the nearest future there is no 
opportunity to observe the H triangles in practice. 


~0(A2), 



The diagrams for the decay 





+ 

C 


are shown in Fig. [3 where the negative relative sign caused by the Pauli interference should 
be taken into account. The corresponding color factors are given by and in ©• 



5 ° 


~‘bc 







Figure 7; The diagrams of 6-quark decay contributing to the weak transition > D^E'^. 


The consideration of decays into the meson is more complicated because of the weak 
scattering of constituent 6 and u quarks as shown in Fig. |H1 We can easily hnd that the 
gluon emission from the other quark lines is suppressed by the color factor since such the 
exchange by the gluon leads to the baryon color-structure factor (Fs in contrast to the (7^ in 
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Figure 8: The diagrams for the weak scattering of b and u quarks contributing to the transition 


~‘bc 


the diagrams shown in Fig. |S1 The spectator charmed quark cannot emit the virtual gluon, 
since the quark line should be on mass shell up to the small virtualities about the relative 
momentum of the quark inside the baryon. 

Further, the color factors for the diagrams in Fig. [Hlhave the form 




rsj 


1 

7K' 


Thus, these factors are of the same order of magnitude as for the decay amplitudes shown in 

Fig. P 





~‘bc 





Figure 9: The diagrams of 6-quark decay contributing to the weak transition ^ 


Finally, in this section we have analyzed the color and weak-interaction structures of decay 
amplitudes and isolate those of the largest magnitude, while an illustrative numerical estimate 
is presented in the next section. 

3 Numerical estimates 

In this section we formulate the framework of a potential model, which allows us to evaluate 
the characteristic widths and branching ratios for the modes under study. 

Let us consider the decay of —>• So, we define the doubly heavy baryon state in 

its rest frame by the following form: 

l-L) = J {vi)C^c{vi)^ d{vi) ■ 

al[b] {k) a] [c](-fc) al[d] (g) |0), 
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( 10 ) 





















where C is the charge-conjugation matrix, and we have introduced the following notations: 


• the relativistic normalization-factor y/2Mi with Mi being the baryon mass in the initial 
state, 

• the wavefunction \I/ given by two factors of the doubly heavy diquark and the d quark, 

• the color wavefunction determined by y/Nj, 

• the spin wavefunction of scalar doubly heavy diquark as determined by the spinor factor 

75/V2, 

• the quark spinors depending on the four-velocity of the baryon, vi, as normalized by 
the condition of sum over the polarization states, say, 

J2Kvi)Kvi) = ^(1 

• denoting the creation operator marked by the flavor, color and momentum in an 
appropriate way. 

In (jl()|l . we use the momentum-space wavefunction in the form of 

kn 


where the four-vector satishes the condition 


= 


= f?i] 




exp 


k ■ vi = 0. 

The model parameter ut is related with the wavefunction T(0) at the origin in the conhguration 
space. Such the quantities were calculated by solving the Schrodinger equation with the static 
potential as described in [T2], so that 

^6c(0) = 0.73 GeV^/^ §^(0) = 0.53 GeV^/^ 

while for the baryon we take 

§,,(0) = 0.61 GeV^/^ ^rf(O) = 0.53 GeV^/^. 


We will see that for the mode under consideration the absolute values of above parameters 
are not so critical, while their ratios are important. 

Then the matrix element is equal to 


Af) — 


Gf 

V2 


V:,V,s ^d(O) y/Ka2 -T-O- u{v2)u{vi), 


( 11 ) 
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where the mass and wavefunction factors as well as the spinors and fonr-velocities are trans¬ 
parently denoted, while T is given by the spin structnre of the matrix element 

^ + ?^ 2 ) 7 /.(l - 75)75(1 + - 75)] , 

and O represents the overlapping of wavefnnctions, which can be calculated in the model 
specihed. So, it is presented by the product of ^ factors 

d = ^bc^cs{y) ■ ^d{y)^ 

where the diquark transition is determined by 


^bc^cs 


‘2uJbc^c 


^^bc^c 


2, ,2 


^hc + ^cs V ^cs + y 


exp 


be 


m?{y‘^ — 1 ) 

L ^1 + y‘^^lc 


( 12 ) 


with y = vi ■ V 2 and 


Mi + rus - mb 
m = me - 

me + mg 


^ me+ 0 {ms^d/me) • 


In the given kinematics the product of four-velocities is hxed by 


+ Ml - Ml 
^ 2 M 1 M 2 

The factor of ^d{y) has the form similar to (j12jl under the appropriate substitutions for the 
values of uj as well as for m —>■ md m- 

Having written down eq. (dD, we have neglected the a^-corrections following from the 
diagrams shown in Fig. El This approximation is theoretically sound. Indeed, the gluon 
virtuality is determined by the expression 


fcj = ml{v 2 + VDY = m\ 


Ml - {M 2 - Mpf 
2,M2 M]j 


{mb + mef + 0{mg^d/me,b) > A: 


QCD- 


The suppression factor of appropriate diagrams is given by 


S ~ 


i^(o)r 


ag{kl) Ne 

kl AEq’ 


where 

|T(0)p ~ Aqcj3 

is the characteristic value of wavefunction, and the virtuality of heavy quark line connected 
to the virtual gluon is of the order of 


AEq ~ me,b- 


10 













Therefore, 


1 


5 


A3 

^'^QCD 


< 1 . 


m3^ In me, fe/AqcD 

Similar relations can be found, if we consider the value of Os-correction due to the weak 
scattering shown in Fig. |S1 One can easily find that the corresponding suppression is given 

by 


5 


Aqcd 

rrich 


< 1 , 


where we have suggested the constituent mass of light quark 

TTT^u,d ~ AqcD- 

Therefore, the above Os-corrections can be neglected to the leading order in the l/mq- 
expansion. 

Further, 

_ , Ml- M 2 

T = 2(1 + 1/) 


M, 


D 


and we can get the matrix element squared 


,2 


\-^d\" = + Ki,yc,\ - Mj)"(l + j,)’ fl-O'y [N^a^r, 


(13) 


where we have expressed the wavefunction of D meson in terms of its effective leptonic con¬ 
stant /d. 


/d — 2 


M, 


»( 0 ), 


D 


Then, the width is given by 


Wl 


-OnO 


D^] = 


\k 




IGttM? 


where fen is the momentum of D meson in the SA rest-frame. 


of) 10 LliC liiUliiClilj Ulii Ui liiCDUli Hi LliC 

Numerically, at Mi ^ 6.9 GeV jinilll] we hnd 


.j, ^ a 1.3 ■ 10-» ps-' X 




ub 




0.0032 0.9752 


L 


D 


(0.222 GeV)2 


■ ■ 


{N,a2f 


Putting nnni 
we get 


^bc. 


= 0.27 ps. 


B[r^ 


be 


: 0 n 0 


D^] ^ 3.6 ■ 10"® X 


VI 


V 


0.0032 0.9752 


/. 


D 


(0.22 GeV)2 


■ C >2 




(14) 


(15) 
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The calculation of overlap between the wavefunctions is mo del-dependent, though one can 
expect that (9 ~ 1. In the framework of potential model described we get 

O ^ 0.44, 


which gives our hnal estimate 

^ ^ 0.7-10"® X 




0.0032 0.9752 


f. 


D 




(0.22GeV)2 0.2 


(16) 


Therefore, we could expect that the other branching ratios are of the same order of magnitude. 


4 Conclusion 

In this work we have extended the reference-triangle ideology to the model-independent ex¬ 
traction of CKM-matrix angle 7 from the set of branching ratios of doubly heavy baryons 
exclusively decaying to the neutral D mesons. Tagging the flavor and CP-eigenstates of such 
the D mesons allows one to avoid the uncertainties caused by the QCD dynamics of quarks. 

We have estimated the characteristic branching ratios in the framework of a potential 
model, which yields, for example, 

^ S“£>»1 « 0.7 ■ 10-». 

Accepting the above value, we can estimate the rate of events at the LHC collider in the 
experiment LHCB or in the BTeV facility at FNAL. So, the production cross section yields 
the characteristic value of about 10® doubly charmed baryons per year [2]. Next, the estimated 
branching rations for the charmed strange baryons are measured for so that the detection 
of charged particles in the final state would cover about 50% of decay events with which 
we accept for the optimistic estimates. The efficiency of observing the neutral charmed meson 
crucially depends on the possibility for detecting the neutral kaons and pions. So, removing the 
neutral kaons and pions could kill the opportunity offered in the present paper. We consider 
the optimistic case with the required detection at work, which gives the efficiency for observing 
the decay at the level of 25%. Therefore, putting the vertex reconstruction efficiency equal 
to 10%, we can expect the observation of about 10® ■ 0.7 ■ 10“® ■ 5 ■ 10“^ ■ 2.5 ■ 10“^ ■ 10“^ ~ 11 
events per year. However, the situation is less optimistic for the detection of CP-eigenstates 
of D^. Indeed, the CP-even mode with the or K~^K~ final states covers only a 5 ■ 10“^ 

fraction of H® decays, which makes the observation unreachable. The CP-odd events with 
Kstt^ can be detected with the branching fraction about 2.5 ■ 10“^, which downs the rate to 
1 event per year. Thus, the reconstruction of reference triangle after 10 years of data taking 
would give, at least, the 30% accuracy for the characteristic triangle side, which makes the 
extraction of angle 7 rather academic exercise in the method offered. 
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